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Abstract 

Solar and atmospheric neutrino anomalies are described by bimaximal mixing of three 
Majorana neutrinos. Neutrino oscillations in appearance fe ^ i^^, Ur and in disappearance 
Vf, long baseline and atmospheric experiments are sensitive to deviations from the ideal 
bimaximal mixing. It is suggested that these deviations may be dominated by a rotation 
in the electron-muon plane of the generation space. Simple seesaw models are described 
supporting this idea. The rotation angle is estimated from Fritzsch's relation. Predictions 
are presented for oscillations in long baseline experiments, for solar neutrinos, and for the 
rates of neutrinoless double beta decays. 
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1 Introduction 



Neutrino masses and oscillations belong to the most exciting area of particle physics. 
It is exactly this area where the Standard Model of quarks and leptons is challenged in a 
most successful way. The recent discovery of muon neutrino oscillations Q and the evidence 
for solar neutrino deficits]^ in comparison to the Standard Solar Model demonstrate that 
the neutrinos are massive. Phenomenological analyses favor two solutions of the solar and 
atmospheric neutrino problems, see e.g.[B|-||T3| and for a recent review see 



large z/^ — Vr mixing for atmospheric anomalies and matter enhanced [|l^ small mixing 
angle oscillations (SMA MSW) for solar neutrinos, see e.g. P^. 



• vacuum oscillations and bimaximal|]l3^|^ or nearly bimaximal||2^ mixing of three light 
neutrinos. 

We believe that one of these two possibilities is true although other options are not com- 
pletely excluded, see e.g. 0, 0, In a near future experimental information from Su- 
perKamiokande and SNO will select the best scenario. In fact a recent measurement of day- 
night asymmetries [p9| ] imposes constraints on SMA MSW solutions. On the other hand the 
vacuum oscillation solutions may be confirmed or excluded by final results on seasonal effects 
and recoil energy distributions, see |]lT], Q and references therein. A preliminary data from 
Super Kamiokande PT|] are analyzed in Ref.[|I^. 

In this article we discuss some consequences of the bimaximal mixing scenario. Evidently 
solutions of 'just so' type[^ are preferred by the existing data. However we do not exclude 
a possibility |]18| that a larger range of Amg is allowed for solar neutrinos if errors are un- 
derestimated in one of the radiochemical measurements. In such a case 'just so' fits are also 



better [11, 13 1. In the following discussion we use as input parameters] 33] 



for oscillations of solar and 



Am^ = 4.3 X 10 



-10 



eV^ 



1.5 X lO"'' eV^ < ^m^atra < 6 X lO"'' eV^ 
for oscillations of atmospheric neutrinos. 

2 Maki— Nakagawa— Sakata lepton mixing matrix 

The Dirac masses of quarks and charged leptons exhibit a strongly ordered hierarchical structure 



m. 



We assume a similar hierarchical structure for the Dirac masses of neutrinos. However this 
hierarchical structure is drastically modified by huge Majorana masses of righthanded neutrinos. 
As a result of the seesaw mechanism||3^ three light and nearly lefthanded neutrinos appear 
in the mass spectrum. Their Majorana masses can form patterns very different from the 
hierarchical one. For the same reason the structures of quark and lepton mixing matrices are 
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also quite different. For quarks the Cabibbo-Kobayashi-Maskawa matrix 1 35] Vckm is nearly 
diagonal and its largest off-diagonal elements 



\Vus\ ^ \Vcd\ ~ ~ 0.22 (1) 

are fairly small. All other off-diagonal elements of Vckm are much smaller and can be parame- 
terized by higher powers of the Cabibbo angle 9c ||3^. Vckm can be written as a product of two 
matrices describing unitary transformations of quarks with weak isospin projections I3 = +1/2 
{u,c,t) and I3 = —1/2 {d,s,b) respectively: 

VcKM = V+^V. . (2) 

In the Standard Model only Vckm is observable. In models of quark masses and mixing both 
matrices V± are specified and related to Yukawa couplings of up and down type quarks. A 
particularly interesting approach [^] relates angles in V+ and V- to mass ratios of up and down 
type quarks respectively. In this way a relation is derived for the mixing between the first two 
generations 

Or 







Ue2 


Ue3 














Ur2 


Ur3 
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— Umns 


V2 











\J nid/ TTis + e^'^^J vfiul vfic ^ \Jmd/ms . (3) 

The mixing angle 9c is dominated by the contribution from the = —1/2 sector. 

The flavor mixing matrixp| for three light neutrinos relates the neutrino mass and flavor 
eigenstates: 

f Uc\ ( Ucl Uc2 Uc3\ f ( Vl\ 

(4) 

Umns can be written as a product of two matrices U± describing transformations of the neu- 
trinos {I3 = +1/2) and the charged leptons {I3 = —1/2): 

Umns = UJU+ . (5) 

It is plausible that the pattern of mixing angles in f/_ describing the charged lepton sector 
resembles that in V- for the down type quarks, i.e. there is an appreciable mixing in the 
electron-muon plane whereas the other mixing angles are much smaller 



f/_ ^ s c (6) 




where [p7 1 



\s\ ~ ^me/m^ ^ 9c/?, ~ 0.07 and c = Vl - ^ 1 - s"^ /2 . (7) 

The second equality in (0) follows from grand uniflcation relations p8 

rud ~ 3me , rris ^ m^/3 . (8) 

In the sector Is = +1/2 the situation is quite different. The structure of the matrix U+ is 
not directly related to the Dirac mass matrix nij-, for neutrinos. It is strongly affected by the 
structure of the Majorana mass matrix Mj^ for the righthanded neutrinos. As a consequence 



2 



the structure of [/+ is different from the structure of 
and atmospheric neutrinos is close to bimaximal, 



We assume that the mixing of solar 



^ U, 



brn 



1/^/2 1/^/2 
-1/2 1/2 l/y/2 
1/2 -1/2 1/^2^ 



(9) 



We also assume that the dominant contribution to deviations from the ideal bimaximal mixing 
originates from the sector = —1/2, i.e. from the form of the matrix [/_. (Later in this paper 
we will demonstrate simple seesaw models in which corrections to eq.(^ are fairly small.) The 
main goal of the present paper is to describe phenomenological consequences of these two 
assumptions. For the sake of simplicity we neglect CP violation and consider all elements of 
Umns as real. Then the lepton mixing matrix is completely specified. From eqs.(||), @ and 
(y) we obtain 

/c/V2-s/2, c/V2 + s/2, s/V2\ 
Umns ~ -c/2 - s/V2 , c/2 - s/V2 , c/V2 (10) 
V 1/2 -1/2 I/V2J 

with s and c given in eq.(|^. 



3 Neutrino oscillations 

It is well known[|, [T^ that for Ami, -C Am^^^ and small {Ug^l, oscillations of solar neu- 
trinos effectively decouple from oscillations in long baseline (LBL) experiments and those of 
atmospheric neutrinos. Both types of oscillations can be described by effective two-neutrino 
oscillations characterized by phase differences 

A. = ^ (11) 

for solar, and 

^atm - 4^ (12) 

for LBL and atmospheric oscillations. In the above equations L is the distance traveled by a 
neutrino of energy E. Probabilities of transitions between different neutrino flavors in atmo- 
spheric and LBL experiments are given by, see and references therein: 

VlBL {Va ^ Vp) = Aafj Sin^ Aatm (« 7^ P) 

VlBL {I'a ^ l^a) = I- Ba siu^ Aatm , (13) 



where (a, j3 = e, ^, r) 



— 4|f/Q,3|^|?7/33|^ , 



Ba = 4|f/,3|'(l-|t/a3|') • (M) 



Then eq.(|T^) implies 



B, = l-c\ B^ = l-s\ Br = l, 



(15) 
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and using we estimate that for appearance type v^, ^ and v^, ^ LBL experiments 

,2 r)„Q ~ n nnc; (16) 



sin^ ~ 0.005 



whereas 



sin^ Id, ^ 0.01 



(17) 



for disappearance of Vf. and i^e neutrinos. These predictions can be tested in future high precision 
experiments. However the numbers are rather small. The estimation ([T6|) is about two times 
smaller than planned sensitivity of MINOS |^9|] and eq.(0) predicts disappearance of electron 
neutrinos at the level 20 times smaller than the present limit from CHOOZ||0[. Still these 
estimations are larger than those following from SMA MSW solutions for which sin^ ^d^, is driven 
to much smaller values by the recent data on day- night asymmetries for solar neutrinos] 
For solar neutrinos, c.f.0] and references therein, 



(z/e ^ Ve) = (l - If/esi') (l - siu^ 2^9^ siu^ Aq) + \U 



where 



For the mixing matrix ([T0| ) 



sin2 2^9e = 4|f/,inf/,2|V [\-\\J, 



^e3| 



Vq = (\- (l - sin^ sin^ Aq) + C»(s^) 



(19) 



(20) 



with 



sin^ Idr,, 



0.99 



/Q ~ X — ,5 ~ U.»» . (21) 

It is evident that the simple picture proposed in this paper can be falsified when sin^ I^q derived 
from solar neutrino data is significantly below 1. In Ref.|0] the best fit to vacuum oscillations 
is 



6.5 X 10"" eVl 



sin^2'(9^ = 0.75 



(22) 



a new 



which seems to be in conflict with the estimation in (21). In Ref.[13|, see also [11 
preliminary data from SuperKamiokande are included into analysis leading to acceptable 
solutions in three regions (called A, C and D) of parameter space 



C 
D 
A 



Am| 

Ami, 
Am| 



4.4 X 10^1° eV^ 

6.4 X 10-1° eV^ 

6.5 X 10-11 eV^ 



sin^ 2^90 
sin^ 2'dQ 
sin^ 2^9^ 



0.93, 
1.00, 
0.70, 



9of -- 

gof 
gof 



14% 



6% 



(23) 



In eq.(PBD the parameter gof (goodness-of-fit) is the probability that a random repeat of the 
given experiment would observe a greater x^, assuming the model is correct. The region A is 
the same as selected by the fits of Ref. [j^. Clearly the analysis of Ref.|[T^ gives results in better 
agreement with eg. (plf) . It is also seen that new data on solar neutrinos will seriously test the 
estimations presented in this article. 
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4 Neutrinoless double beta decay 

Let us discuss now predictions for neutrinoless double beta decays Probabilities of 

such transitions depend on a mass parameter B which is equal to the absolute value of the 
element (A^a)ii of the matrix 

Na = Umns Ml,a Umns^ . (24) 

Ml^a (for A = 1 1 11 1 1ll, IV) denotes a diagonal matrix of Majorana masses for the three light 
neutrinos: 

Ml^a = m diag{Xi,X2,X3) , (25) 

where |Aj| < 1 + 0(e, ^); e and ^ are small parameters to be discussed in the following section. 
In our convention the eigenvalues Aj are ordered in such a way that 

Aml/m^ = \\l- Xl\ (26) 

and 

Aml,Jm'^\\l-\l\^\\l-Xl\ . (27) 

For bimaximal mixing i/g is an equal mixture of the mass eigenstates 1 and 2. 

Considering Amg/Am^^^ as a small perturbation one obtains |Ai| = |A2|. Then, up to 
irrelevant sign redefinitions there are only four mass patterns consistent with the data on Qv2(3: 



I) m = JAmlt^, Ai = A2 = 0, A3 



II) m = jAmlt^ , Ai = -1, A2 = 1, Ag 



III) m = JA'mlt^, Ai = A2 = 1, A3 



IV) m = C(l eV) > ^Amlt^ , Ai = -1, A2 = 1, A3 = r/ where i] = ±1 . 

In case IV Am^j^/m^ is considered as a small perturbation^. The states 1 and 2 have opposite 
CP parities (i.e. Ai and A2 have opposite signs, see e.g.[|l5) due to the experimental bound[^]: 

B < 0.2 eV (28) 

which is violated for m > 0.2 eV and Ai = A2 = 1. The argument is analogous to the case /// 
which is considered in the following, c.f. eq.(^). 

If corrections to bimaximal mixing are neglected the following mass matrices are obtained 



Nj = WAm2,„ Oil ' (29) 




model independent analysis is presented in a recent preprint|41 . 
^The case Am^^^/m^ = C(l) requires much more tuning of parameters than the case IV. We do not discuss 
this case in the present article because it is not clear if a reasonable model of this kind exists. Furthermore an 
original motivation for degenerate neutrino masses is an appreciable neutrino contribution to dark matter in 
the Universe whereas this contribution is small if Am^^^/m^ is not small. 
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II 




III 



Niv 



where for A = I, II, III, IV: 

In the following section we consider also matrices 

Na = U+ Ml,a 



(30) 
(31) 
(32) 

(33) 
(34) 



which include deviations from the ideal bimaximal mixing due to the sector J3 = +1/2. The 
textures (p^ ) -(|32| ) have been discussed in the literature [|19|-|^. It is seen that only Nm 
leads to non-zero rates of Oz/2/3 decays. In this case deviations from bimaximal mixing can be 
neglected and one obtains the following prediction 



|iV,,,|n ^ V^^'*™ = (6 ± 2) X 10^^ eV . (35) 
In the three other cases the role of [/_ is essential. Keeping only leading terms one obtains 



\Ni 


11 ~ 






= (1.5 ±0.5) X 


\Nn 


11 ~ 






= (6±2) X 10" 


\Niv 


11 ~ 


72 


\s\ m 


^ 0.1 m . 



eV 



(36) 
(37) 
(38) 



Evidently the cases /// and IV can be confirmed or ruled out by next generation experiments 
which may be sensitive to B as low as 0.01 eV[^. Even the present limit implies that for 
degenerate masses m must be smaller than 2 eV. 



5 Patterns of neutrino masses 

In this section effects are considered due to non-zero value of the ratio AmQ/Am^^^ which 
were neglected in the discussion of the preceding section. We present simple seesaw models 
which show that these corrections are quite small, so the estimations (P^-(P5D are not much 
affected. In the following we consider the four cases I-IV corresponding to different neutrino 
mass patterns. 



5.1 Case I 



A class of seesaw models corresponding to case / has been described in|jT9|. In these models the 
third mass eigenstate is much heavier than the other two whose masses are fairly close. In the 
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following this mass pattern (mi fsi m2 -C ^3) is called semi-hierarchical. In terms of a small 
parameter 



Ami V^' 



0.004 (39) 



the eigenvalues Aj are given by the expansions: 

Ai = e-eV2 + ... 
A2 = -e-eV2 + ... 

A3 = l + e2 + ... (40) 

For semi-hierarchical mass pattern z/g is a linear combination of the two lighter mass eigenstates 
of opposite CP parities. The element {Nj)ii of the matrix 

Ni = U+ Ml,i f/+^ (41) 

is strongly suppressed and the main contribution to 01/2/3 transitions originates from the mixing 
in the sector J3 = —1/2. The estimation of lA^^/ln given in eg. (^6]) is not affected as can be seen 
from the explicit form of the matrix f/+ :0 

/(l + e/2)/v^ (l-e/2)/v^ e \ 
f/+ = -1/2 - e 1/2 - e 1/v^ + 0{e') = U^m + 0{e) (42) 

V 1/2 -1/2 I/V2J 

and 

_ / 2^/\\ 

Nj = l\ 1 1 + 0(e2) . (43) 
\2'/^e 1 1 J 

In case / a strongly ordered hierarchical mass pattern (mi <^ m2 <^ m3) is obtained if e.g. 

Ai = + . . . , A2 = e + . . . , A2 = 1 + . . . , (44) 

with 

|e| = ^Aml/Aml,^ ^ 3.7 x 10"^ . (45) 

The magnitude of e is an order of magnitude smaller than for semi-hierarchical spectrum. There 
is no suppression of (^^7)11 because the lighter states have the same CP parity but the parameter 
e is small. Seesaw models leading to hierarchical spectrum may require some fine tuning. Thus 
we believe that hierarchical mass pattern is less attractive from theoretical point of view. In 
spite of theoretical prejudices we note that for both patterns predictions for oscillations and 
Oz/2/3 decays are similar. 



In notation of Ref. |19[| r — l?l'^t. The columns of the matrix J7+ in eq.(|4^) are normahzed eigenvectors ui, 



W2 and W3, c.f. eq.(29) in |19 . We have corrected a misprint in the second element of the vector V2 in [|H 
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5.2 Case II 



Corresponding seesaw models can be easily derived. An example is (in notation of UToi 



m ~ 

D 



/x^y 
X —X 
V 1 



and M"^ 



^ ail CLi2 
ai2 

V 0, 



(46) 



where nij-, is the Dirac mass matrix and M denotes inverse of the Majorana mass matrix for 
the righthanded neutrinos^. Then 



Nji = U+ Ml J J = mjm-^^ = 



atm 



2e 1/v^ 
1/V2 
-1/V2 



(47) 



with the eigenvalues 



Ai 



Vl + e2 + e 



A, 



Vl + e2 + e 



and 



|e| = 2-3/2^y|an/ai2| = Amy(4Am^,^) = 3.4 x 10 



(4J 



(49) 



Evidently corrections of order e to Umns are negligibly small. The spectrum of light neutrinos 
contains two heavier states. The electron neutrino is a nearly maximal mixture of these states 
with a tiny admixture of the third state which is massless. We call this spectrum semi- degenerate 
(mi 7712 ^ ^s)- The ratio |aii/ai2| must be very small which implies that the two heavy 



Majorana neutrinos form a pseudo-Dirac systemj^. In our opinion the semi-degenerate mass 
pattern is more attractive for solutions of the solar neutrino problem with larger values of 
Amg/Am^j^ [|l^, see also discussion in p4|. 



5.3 Case III 

As in case II the parameter e is tiny and the related corrections are negligible. 



5.4 Case IV 



For degenerate neutrino masses 1 22, 23] one can choose the eigenvalues Aj in the following way: 



Ai 



-1 + e 



with 



|e| = AmQ/(4m^) and 
Corrections to eg . (^8)) are small for |^| ^ 



A2 = l + e, A3 = ?7 + ^, 
lel = A7nlJi2m') 



(50) 
(51) 



This form of M^-"- means that either one of the righthanded neutrinos is much heavier and decoupled or 
there are only two heavy neutrinos in the particle spectrum and their couplings to the lefthanded neutrinos are 



given by the matrix 



xy 
1 -1 
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If the case of degenerate neutrino masses is a reasonable option for vacuum oscillations 
remains an open question. On one side dynamical models have been recently found |^6|] for the 
degenerate mass pattern. On the other side the parameter e is extremely small and the problem 
of stability may be very serious, see It is also not clear if neutrino masses in the range of 
1 eV are really needed by cosmology, see e.g. Fortunately it is exactly this mass range and 
mass pattern for which our estimations give numbers close to the present experimental limit 
for decays. Therefore one may hope that a definite answer will come from the experiment 
in not very distant future. 



6 Summary 

Neutrino oscillations and Majorana masses are discussed assuming bimaximal mixing of leptons. 
In this scheme the atmospheric neutrino anomaly is described as a result of Vj. 
oscillations and deficits of solar neutrinos [^] as vacuum oscillations of i^e- A simple hypothesis 
is proposed that deviations from bimaximal mixing are dominated by a rotation in the electron- 
muon plane of the charged lepton sector = —1/2. This hypothesis is in agreement with the 
present data and can be tested by future precision data on solar neutrinos. Predictions for 
oscillations in long baseline experiments are given. There are four classes of neutrino mass 
patterns leading to quite different predictions for the rates of neutrinoless double beta decays. 
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